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We consider tlie classical and quantum dynamics of excitations in a system of two capacitively 
coupled Josephson junctions. In the classical case the equations of motion admit discrete breather 
solutions, which are time periodic and localized predominantly on one of the junctions. In the 
quantum case breather states are found in the central part of the energy spectrum of the confined 
nonescaping states of the system. We perform a systematic analysis of their tunneling frequency, 
site correlations, fluctuations of the number of quanta, and entanglement. Quantum breather states 
show strong site correlation of quanta and are characterized by a strong excitation of quanta on 
one junction which perform slow coherent tunneling motion from one junction to the other. They 
suppress fluctuations of the total number of excited quanta. Quantum breather states are the least 
entangled states among the group of eigenstates in the same range of the energy spectrum. We 
describe how quantum breather excitations could be experimentally observed by employing the 
already developed techniques for quantum information processing using Josephson junctions. 



PACS numbers: 63.20.Pw, 74.50.+r, 85.25.Cp, 63.20.Ry 



I. INTRODUCTION 

Josephson junctions are the subject of extensive stud- 
ies in quantum information experiments because they 
possess two attractive properties: In their classical 
regime they are nonlinear devices, but also show macro- 
scopic quantum behavior [H, d, 0] . The dynamics of a bi- 
ased Josephson junction (JJ) is analogous to the dynam- 
ics of a particle with a mass proportional to the junction 
capacitance Cj, moving on a tilted washboard potential 



U{ip) = -Ic— cos ip~ hip— , 



(1) 



which is sketched in Fig[T]-b. Here ip is the phase dif- 
ference between the macroscopic wave functions in both 
superconducting electrodes of the junction, If, is the bias 
current, Ic is the critical current of the junction, and 
$0 = h/2e the flux quantum. When the energy of the 
particle is large enough to overcome the barrier AC/ (that 
depends on the bias current /;,) it escapes and moves 
down the potential, switching the junction into a resis- 
tive state with a nonzero voltage proportional to ip across 
it. Quantization of the system leads to discrete energy 
levels inside the potential wells, which are nonequidis- 
tant because of the anharmonicity. Note that even if 
there is not enough energy to classically overcome the 
barrier, the particle may perform a quantum escape and 
tunnel outside the well, thus switching the junction into 
the resistive state [l|. Thus each state inside the well is 
characterized by a bias and state-dependent lifetime, or 
its inverse — the escape rate. 

Progress on manipulation of quantum JJs includes 
spectroscopic analysis, better isolation schemes, and si- 
multaneous measurement techniques 0, H, 0, H, [E 13' 
and paves the way for using them as Josephson-junction 
qubits in arrays for experiments on processing quantum 
information. Typically the first two or three quantum 



levels of one junction are used as quantum bits. Since 
the levels are nonequidistant, they can be separately ex- 
cited by applying microwave pulses. 

However, improvements in experiments manipulating 
Josephson-junction qubits may have applications beyond 
the processing of quantum information. Operating the 
junctions at larger energies in the quantum regime may 
give rise to other interesting phenomena that nowadays 
can be experimentally observed by using already devel- 
oped experimental techniques. For instance, it was sug- 
gested that JJs operating at high energies may be used 
for experiments on quantum chaos [ol, [loj . 

Another interesting phenomenon is the excitation of 
discrete breathers. They are time periodic space local- 
ized excitations in anharmonic lattices with translational 
mvariance [O, [H, [H, El- They locahze energy expo- 
nentially in space for short-range coupling between lat- 
tice sites, and have been experimentally observed in such 
different systems as bond excitations in molecules, lat- 
tice vibrations and spin excitations in solids, electronic 
currents in coupled JJs, light propagation in interact- 
ing optical waveguides, cantilever vibrations in microme- 
chanical arrays, cold atom dynamics in Bose-Einstein 
condensates loaded on optical lattices, among others 

[H, HE [13, H [H, H mi [13, m . 

In the quan tum regime, quantum breathers (QB) 

MMJ^iM m, m, m, IsH, m m, m, HE HE [33, 

ISSl . [39l . HE EH appear as nearly degenerate many-quanta 
bound states. Though being extended in a translation- 
ally invariant system, they are characterized by expo- 
nentially localized correlation functions in full analogy to 
their classical counterparts [lEIH]- When such states su- 
perpose the result is a spatially localized excitation with 
a very long time to tunnel from one lattice site to an- 
other. At variance to the classical case, the evolution 
of these excitations in time has not been experimentally 
studied in detail. So far they have been indirectly ob- 
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served 



spectroscopic analysis in molecules and solids 

. I, In, 113,141, mi, Hi- 

The possibility to directly observe QB excitations 
evolving in time was addressed by us in a letter [sl] for a 
system of two capacitively coupled JJs, where by calcu- 
lating the eigenstates and the spectrum of the system we 
identified QB states as weakly splitted tunneling pairs 
of states [31I, m, m, [131 ■ These eigenstates appear in 
the middle of the energy spectrum of the system and 
are characterized by correlations between the two junc- 
tions - if one of them is strongly excited, the other one is 
not, and vice versa. By exciting one of the junctions to 
a large energy (many quanta), we strongly overlap with 
QB tunneling states. Consequently we may trap the ex- 
citation on the initially excited junction on a time scale 
which sensitively depends on the amount of energy ex- 
cited, and on the applied bias. We described how QB 
excitations could be directly observed in time using the 
available techniques for manipulating JJs in the quantum 
regime. 



In this work we present an extended analysis of the 
system, performing a systematic and comparative anal- 
ysis of different properties of QB states. We study their 
tunneling rates, the site correlations of excited quanta, 
the fluctuation of the number of excited quanta, and the 
entanglement of the QB states. 

In section II we describe the model for the two coupled 
JJs [13, [1^ [1^ and briefly consider the classical dynam- 
ics, where the equations of motion are numerically solved 
finding discrete breather solutions. In section III wc con- 
sider the quantum model and introduce the basis wc use 
to numerically diagonalize the Hamiltonian matrix. We 
define correlations functions which, together with the en- 
ergy spectrum, will help us to identify QB states. Then 
we compute the time evolution of initially localized ex- 
citations and relate it to the spectral properties of the 
system. In section IV we address the fluctuation of the 
total number of quanta in the eigenstates. In section V 
we explore the entanglement of the eigenstates. In sec- 
tion VI we describe how QB excitations evolving in time 
could be experimentally observed, and discuss how es- 
caping and decoherence (effects that are not taken into 
account in the quantum model) would affect the obser- 
vations. We conclude in section VII. 



II. THE MODEL AND CLASSICAL DYNAMICS 

The system is sketched in FigUJa: two JJs are coupled 
by a capacitance Cc, and they are biased by the same 
current /{,. The strength of the coupling due to the ca- 
pacitor is C = Cc/{Cc + Cj). The Hamiltonian of the 
system is 



H 



Pi 
2m 



p2 

2m 



C 




(c) 



FIG. 1; (a) Sketch of the two capacitively coupled Josephson 
junctions, (b) Sketch of the washboard potential for a single 
current-biased Josephson-junction. (c) Circuit diagram for 
two ideal capacitively coupled Josephson junctions. 



where 



m = Cj(l + C) 



Pi. 



{Cc + Cj) 



£0 

27r 

£0 
27r 



(3) 
(4) 



Note that the conjugate momenta Pi^2 are proportional 
to the charge at the nodes of the circuit (which are la- 
beled in FigHJ-c). When the junctions are in the su- 
perconducting state, they behave like two coupled an- 
harmonic oscillators with plasma frequency Wp(7) = 

V27r/e/$oCj(l-HC)[l - 7']'/^ 7 = h/Ic being the nor- 
malized bias current. The classical equations of motion 
are given by 



$0 'I'o 
(sin(/3i,2 -I-Csin(p2,i) + 7^ (1 + C)7 • (5) 



m 



P1P2. (2) 



Despite being invariant under permutation of the junc- 
tion labels, these equations admit discrete breather so- 
lutions [l^l, which are time periodic and for which the 
energy is localized predominantly on one of the junctions 
(Figl2). These orbits can be numerically cornputed with 
high accuracy using Newton algorithms [13, [111 . 

The existence of discrete breathers is possible because 
the anharmonicity in the JJ potentials makes the fre- 
quency of these excitations (and all of their harmonics) 
nonresonant with the normal modes — \/l ± C,ujp{'^) 
of the coupled-junctions system, whose corresponding or- 
bits are delocalized [l3|. For the parameters 7 = 0.99 
and C = 0.1, the normal- mode frequencies are w+ = 
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FIG. 2: (a) Time evolution of the phase differences, and (b) 
corresponding sketch, of an exact discrete breather solution of 
the equations of motion ([5]) with frequency ili, — 0.303ajp(0) 
(time is measured in units of the inverse plasma frequency at 
zero bias ujp{0)~^). The parameters are 7 — 0.99 and ^ = 0.1. 



0.394a;p(0) (in-phase mode) and uj- = 0.356tJp(0) (out- 
of-phase mode). The periodic solution shown in Fig[2] 
has a frequency below the out-of-phase mode frequency, 
thus the discrete breather solution is out of phase as well. 



III. QUANTUM DYNAMICS: EXCITING 
QUANTUM BREATHER STATES 

In the quantum case we compute the energy eigen- 
values and the eigenstates of the system. Since we are 
interested only in the energy transfer between the junc- 
tions, we neglect quantum escape for states which will 
not escape in the classical limit. Thus we use a changed 
potential energy for the single JJ by adding a hard wall 
which prevents escape: 



U {(p) if (fi < TT — ipo 
00 if if > TT — ifQ 



(6) 



where ipo = arcsin7 is the position of the minimum of 
the potential and ir — ipo gives the position of the first 
maximum to the right from the equilibrium position ipo 
(FigdJ-b). We will later compare the obtained tunneling 
times with the true state dependent escape times. 

The Hamiltonian of the two-junctions system is given 

by 



H = Hi+H2 + CV, 



(7) 



where Hi ~ /2m + Uq{(pi) is the single-junction Hamil- 
tonian and V = PiP2/m is the interaction that couples 
the junctions. The eigenvalues Sm and eigenstates \ni) 
of the single-junction Hamiltonian Hi were c omp uted by 
using the Fourier grid Hamiltonian method [59(. \ni) is 
also an eigenstate of the number operator fii with eigen- 
value Hi. In the harmonic approximation f6d\ 



(8) 



where aj and a-i are the bosonic creation and annihilation 
operators. Since only states with energies below the clas- 
sical escape energy (barrier) are taken into account, the 



computed spectra have a finite upper bound. The pertur- 
bation V does not conserve the total number of quanta 
ni + n2, as seen from the dependence of the momentum 
operators on the bosonic creation and annihilation oper- 
ators in the harmonic approximation: 



Pi, 



{^0/2TT)^{l+C)CjhLOp/2 



X (a 



1.2 



(9) 



The Hamiltonian matrix is written in the basis of prod- 
uct states of the single-junction problem {\ni,n2) = 
\ni) (8) |?^2)}. The invariance of the Hamiltonian under 
permutation of the junction labels allows us to use sym- 
metric and antisymmetric basis states 



\ni,n2)s,A 



1 



V2 



{\ni,n2) ± 1^2,^1)) 



(10) 



to reduce the full Hamiltonian matrix to two smaller sym- 
metric and antisymmetric decompositions of H, which af- 
ter diagonalization respectively give the symmetric and 
antisymmetric eigenstates of the system. 

In order to identify quantum breather states, whose 
corresponding classical orbits are characterized by energy 
localization, we define the correlation functions: 



/^(l,2) = (nin2), 



(11) 



(12) 



where (A)^ = (XmI^IXm>> {IX/^)} being the set of eigen- 
states of the system. The ratio < /^(l, 2)//^(l, 1) < 1 
measures the site correlation of quanta: it is small when 
quanta are site-correlated (i.e. when many quanta are lo- 
cated on one junction there are almost none on the other 
one) and close to one otherwise. 

In FiglS] we show the nearest neighbor energy spac- 
ing (tunneling splitting) and the correlation function of 
the eigenstates. For this, and all the rest, we used 
Ic = 13.3 nA, Cj = 4.3 pF, and C = 0.1, which are 
typical values in experiments. We see that in the cen- 
tral part of the spectrum the energy splitting becomes 
small in comparison to the average. The corresponding 
pairs of eigenstates, which are tunneling pairs, are site 
correlated, and thus QBs. In these states many quanta 
are localized on one junction and the tunneling time of 
such an excitation from one junction to the other (given 
by the inverse energy splitting between the eigenstates 
of the pair) can be exponentially large and depend sen- 
sitively on the number of quanta excited. 

Note that the tunneling of quanta between the JJs 
occurs without an obvious potential energy barrier be- 
ing present (the interaction between the junctions is 
only through their momenta, as seen in the Hamiltonian 
([7])). This process has been coined dynamical tunnel- 
ing [61I [6^ . [63j , to distinguish from the usual tunneling 
through a potential barrier. In dynamical tunneling, the 
barrier — a so-called invariant separatrix manifold — is 
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FIG. 3: (a) Energy splitting and (b) correlation function vs. 
energy of the eigenstates of the two-junctions system (open 
circles, symmetric eigenstates; filled circles, antisymmetric 
eigenstates). The labeled arrows mark the energy correspond- 
ing to the peak of the spectral intensity in Fig|4]-b, d, and f 
(see text). The parameters are 7 — 0.945 and — 0.1 (22 
levels per junction) . 



only visible in phase space, where it separates two re- 
gions of regular classical motion between which the tun- 
neling process takes place. Therefore, when referring to 
the tunneling between the JJs, we implicitly mean that 
it is dynamical. 

The fact that the strongest site correlated eigenstates 
occur in the central part of the energy spectrum may be 
easily explained as follows: Let N be the highest excited 
state in a single junction, with a corresponding maximum 
energy AC/ (FiglT|). For two junctions the energy of the 
system with both junctions in the iV-th state is 2AU, 
which roughly is the width of the full spectrum. Thus 
states of the form |7V, 0) and |0, N) that have energy AU 
are located approximately in the middle. 

Having the eigenvalues and eigenstates, we compute 
the time evolution of different initially localized excita- 
tions, and the expectation value of the number of quanta 
at each junction {hi){t) = {t)\fii\^ [t)) . Results are 
shown in Fig|3^, c, and e. We also compute the spectral 
intensity — \{x^J.\^ o)\^ ^ which measures the strength 
of overlap of the initial state |5'o) with the eigenstates. 
Results are shown in Fig|4]-b, d, and f, where we can see 
a peak in each case, which corresponds to the arrows in 
FigOl-b. We can see that the initial state |^'o) = |0, 5) 
overlaps with eigenstates with an energy splitting be- 
tween them being relatively large and hence the tunnel- 
ing time of the initially localized excitation is short. For 
the case |^'o) = |0,19) QBs are excited: The excitation 
overlaps strongly with tunneling pairs of eigenstates in 
the central part of the spectrum, which are site corre- 
lated and nearly degenerate. The tunneling time of such 
an excitation is very long, and thus keeps the quanta lo- 
calized on their initial excitation site for corresponding 
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FIG. 4: Time evolution of expectation values of the number of 
quanta at each junction (left panels) for different initial exci- 
tations with corresponding spectral intensities (right panels), 
(a) and (b): i*o> = |0,5); (c) and (d): |*o) = |0,19); (e) 
and (f): |vl/o) = |9, 19). Open circles, symmetric eigenstates; 
filled circles, antisymmetric eigenstates. The energies of the 
peaks in the spectral intensity are marked by labeled arrows in 
FigOJb (see text). The parameters are 7 = 0.945 and = 0.1 
(22 levels per junction). 



times. Finally the initial state |^'o) = |9, 19) overlaps 
with weakly site correlated eigenstates with large energy 
splitting. Hence the tunneling time is short again. 

We computed also the time evolution of the expecta- 
tion values of the number of quanta for initial conditions 
which are coherent or incoherent (mixtures) superposi- 
tions of product basis states with equal weights. This 
is relevant for experiments, since it may be hard to ex- 
cite one junction to a determined state but easier to ex- 
cite several states of the junction at the same time. We 
used coherent superpositions (characterized by well de- 
fined states l^'o)), and mixtures (characterized by their 
corresponding density operators po), of four basis states 
around the already used initial states: For the state |0, 5) 
we superposed the basis states |0, 5), |0,6), |0,7), and 
|0,8), for |0,19) the basis states |0,20), |0,19), |0,18), 
and |0, 17), and for |9, 19) the basis states |9, 20), |9, 19), 
|9, 18), and |9, 17). Both for superposition and mixture of 
basis states, the results are qualitatively similar to those 
shown in FigH) Therefore we expect that some impre- 
cision in exciting an initial state in the junctions would 
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FIG. 5: Contour plot of the logarithm of the density of the 

asymmetric state |x) = dxi^') + \Xb^^))/'^ ^ ^ function of 
the number of quanta at junctions 1 and 2 (see text). The 
parameters are 7 = 0.945 and = 0.1 (22 levels per junction) . 



not affect in a relevant way the results, we may also con- 
clude, that the excitation of QB states does not rely on 
the phase coherence. That conclusion will be supported 
later by the study of entanglement. 

Let us estimate how many quanta should be excited in 
the junctions in order to obtain QBs (tunneling pairs). 
We compute the density p(ri 1,712) — |(rii,n2|x)P of 
the asymmetric state \x) = Hxi^^) + \Xb^''))/V^j where 
[x'f '^'') are the eigenstates belonging to a tunneling pair 
|37| . In Fig. [5] we show a contour plot of the logarithm 
of the density for the tunneling pair with energy marked 
by the arrow labeled by number two in Fig. [5]-(b). We 
see that the density has its maximum around ni = 19 
and n2 — which is consistent with the result shown 
in Fig. [4]-c and d where QBs were excited by using this 
combination of quanta in the junctions. 



IV. FLUCTUATIONS OF THE TOTAL 
NUMBER OF QUANTA 

Even though the Hamiltonian does not commute with 
the total number of quanta iV = rii + n2, in Fig[3]-a, c 
and e we see that its expectation value has very small 
fluctuations (less than one). We can see this approxi- 
mate conservation of the number of quanta also in the 
density plotted in Fig|5l where shows a rim along the 
line ni -\- n2 = N (=19). This might be computationally 
advantageous when considering larger systems because 
the strict conservation of = ni -I- n2 would allow us to 
truncate the Hilbert space and work within a subspace 
with fixed N. Each time the interaction operator V acts 
on a basis state with given N, it will generate also states 
with iV ± 2, as can be seen from first order perturbation 
theory in ^ at the harmonic approximation. To study 
these fluctuations numerically we computed the follow- 
ing weight function for each eigenstate, which measures 




N 

FIG. 6: From Left to right: Weight function as a function 
of the total number of quanta N = ni + n2 for the three 
symmetric eigenstates at the peaks of the spectral intensities 
shown in Fig|4]-b, d, and f respectively. The parameters are 
7 = 0.945 and ^ ~ 0.1 (22 levels per junction). 

the relative contribution of all basis states with a given 
N to the eigenstate under consideration: 

W^.{N)^ J2 \{ni,n2\x^)\' (13) 

m ,n2 
ni+n2=N 

In Fig[6]we show the weight function for the three sym- 
metric eigenstates which correspond to the peaks of the 
spectral intensities shown in Fig. \^h, d, and f respec- 
tively. For the lowest-energy state we can see the ex- 
pected appearance of two satellite peaks separated by 
two quanta from the central one. For the higher-energy 
eigenstates the harmonic approximation does not hold. 
Most importantly we see, that for states in the lower and 
middle part of the energy spectrum, the fluctuation of 
the number of quanta is weak, and corresponding states 
contribute less than one percent to the eigenstate. This 
is apparently not true at the upper end of the energy 
spectrum. 

Note that the obtained amplitude of fluctuations in 
the time evolution is much less due to averaging effects 
and the smallness of the strength of the perturbation (. 
The calculation of (N) (t) from perturbation theory in the 
harmonic approximation shows that this quantity oscil- 
lates in time with an amplitude that is proportional to 
C^nomo, where no, and toq are the energy levels initially 
excited in the junctions. For C = 0.1 one finds that the 
fluctuations are of the order of 10^^ in the case shown in 
Fig|3]-a; and of the order of 10~^ in the cases in FigU-c 
and e. Numerical results are consistent with these esti- 
mates. 

To characterize the variation in the total number of 
quanta in the eigenstates we computed the fluctuation: 



^{AN^), = {N)l (14) 
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FIG. 7: Fluctuation of the total number of quanta for the 
eigenstates of the coupled-junction system (open circles, sym- 
metric eigenstates; filled circles, antisymmetric eigenstates). 
The parameters are 7 = 0.945 and C, = 0.1 (22 levels per 
junction). 



In FigHwe plot the relative fluctuation ^(AiV2)^/(iV)^ 
for the eigenstates, where we can see that it is very small, 
and for the QB states in the central part of the spectrum 
it has the smallest values. This follows from the fact that 
QBs are close to eigenstates having the form 



Ix)qb 



V2 



(|n,0)±|0,n)). 



(15) 



with n < N {N is the number of levels per junctions). 
These are eigenstates of the total number operator N — 
fii + ^2, for which the corresponding fluctuations given 
by HI]) vanish. 



V. ENTANGLEMENT OF QB STATES 

Since QB states are close to eigenstates of the C = 
case given by p5|) . one expects that the degree of en- 
tanglement in QB states is similar to the degree of en- 
tanglement in such states. Since only two basis states 
are involved, it can not be a state of maximum entangle- 
ment. We measured the degree of entanglement in the 
eigenstates of the system by minimizing the distance of a 
given state to the space of product states, which depends 
on the larg est eigenvalue of the reduced density matrix 



N 

E 



VAni , 



fni9n2 ) 



(16) 



where Xni.na = ("17"2|Xa')' ^"^"^ functions and 
g„2 are such that A is minimum (see appendix for explicit 
formulas). A measures how far a given eigenstate of the 
two-junction system is from being a product of single- 
junction states, and has values < A < 1 (see appendix). 
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FIG. 8: Entanglement of the eigenstates of the coupled- 
junction system for different values of the coupling strength 
(open circles, symmetric eigenstates; filled circles, antisym- 
metric eigenstates). Here 7 = 0.945. 



For C, — Q the eigenstates of the system are the basis 
states given by eq. (fTO|) . where for rii — n2 it follows 
that A = 0, and for m ^ n2 (which includes the state 
in eq. (|15p ) A = 0.5. This measure has a direct relation 
to the distance of a given eigenstate from a possible one 
obtained after performing a Hartree approximation [6^ . 

In a quantum integrable model with two degrees of 
freedom [3l| it was shown that the region of existence 
of QB states in the energy spectrum is separated from 
other states by the energy threshold for which discrete 
breathers exist in the corresponding classical model. 
Pairs of eigenstates with energies beyond this thresh- 
old show exponentially decreasing energy splitting. In 
a similar quantum model (gtI . [68| , it was shown that at 
the mentioned energy threshold the entanglement (using 
the von Neumann entropy) becomes maximum and then 
decreases with energy. From these two results we expect 
that QB states in our case show decreasing entanglement 
A with energy, tending to 0.5. 

In FiglS] we show the entanglement of the eigenstates 
for different values of the coupling strength For C = 
the entanglement has the values and 0.5 correspond- 
ing to the basis states with equal and distinct number 
of quanta at each junction respectively. When C > 
the eigenstates become linear superpositions of the basis 
states and the entanglement rises, being larger as long as 
more basis states are involved in building up an eigen- 
state. This can be seen in FiglHl where we plot contours 
of the density of four symmetric eigenstates marked by 
labeled arrows in Figl5]for ( = 0.1: the low-energy eigen- 
state marked by the arrow 1 in Fig[8] consists mainly 
of a superposition of a few basis states |ni,n2) fulfilling 
JT-i + "-2 = 3 as seen in FiglS]-a, hence the entanglement 
is relatively small. When going up in energy the entan- 
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(a) (b) 




FIG. 9; Contour plots of the logarithm of the density of the 
symmetric eigenstates marked by labeled arrows in Fig[8]for 
the case C = 0.1: (a) S-5 (arrow 1), (b) S-117 (arrow 2), (c) 
S-95 (arrow 3), (d) S-246 (arrow 4). The normalized bias 
current is 7 = 0.945 (22 levels per junction). 



glement in the eigenstates quickly increases, becoming 
maximum in the central part of the energy spectrum and 
then decreases. An eigenstate in this region of the spec- 
trum hke the one marked by the arrow 2 in Fig [8] involves 
many basis states fulfilling 711+712 = 20 (Fig[5]-b), hence 
the entanglement is large. However, for QB states living 
in the same energy region like the one marked by the ar- 
row 3 in FiglU which has the form shown in eq. (jisp as 
visible in Figl5]-c, the entanglement is smaller and tends 
to 0.5 as expected. Finally, high-energy eigenstates like 
the one marked by the arrow 4 in FiglS] involve not so 
many basis states (Figl9]-d), therefore the entanglement 
is also smaller. 

We also computed the von Neumann entropy [g^ (see 
appendix for explicit formulas), which is another stan- 
dard measure of entanglement, and the results were con- 
sistent with those discussed above. 

We thus conclude, that when QBs appear in a certain 
part of the energy spectrum, their entanglement drops as 
compared to the typical entanglement of nearby states. 
The reason is, that QB states predominantly excite two 
symmetry-related basis states, as opposed to the typical 
excitation of many other basis states. 



VI. POSSIBLE EXPERIMENTAL 
OBSERVATION OF QUANTUM BREATHERS 

The experimental observation of QBs may be possi- 
ble by using the scheme of McDermott et al for simul- 
taneous state measurement of coupled Josephson phase 
qubits Q, where by applying current pulses in the bias 
current through each junction the time evolution of the 



occupation probabilities in the qubits is measured. By 
applying a microwave pulse on one of the junctions we 
excite it into a high energy single-junction state with en- 
ergy £; and leave the other one in the ground state. In 
this way we have an initial state similar to the ones shown 
in Figl4]-c and d. After a variable period of time we ap- 
ply simultaneous current pulses to the junctions to lower 
their energy barriers At/ and enhance the probability of 
tunneling outside the potential well. Then we test which 
junction switches to the resistive state (detected by a 
measurable voltage across it). By repeating the measur- 
ing many times we obtain the populations in the junc- 
tions as a function of the time between the initial pulse 
and the simultaneous measuring pulses. 

Let us discuss the so far neglected quantum escape. 
For that we computed T^scape by using the semiclassical 
formula [70l | 

're~iape(e) = ^^exp|-^y p(^)dV3|, (17) 

where a and h are the turning points of the classical 
motion in the reversed potential at U{lp) — e, p{'^) — 
^j2\U{ip) — e], and a;(e)/27r is the frequency of the os- 
cillations inside the initial well. In table U we show 
the escape time from different metastable states, and 
we compare it with the tunneling time Ttunnel of an ini- 
tial excitation |^(0)) = |0,/) between the two junctions, 
estimated from the energy splitting of the (symmetric- 
antisymmetric) pair of eigenstates with the largest over- 
lap with the initial excitation. We see that for / = 19, 
where we excite QBs, the escape time is long enough 
for observing at least one tunneling exchange between 
the two junctions before escaping to the resistive state. 
Note that the cases / = 18 and 17 also excite QBs which 
would show even more tunneling exchanges before escap- 
ing. The case I = 16 does not excite QBs but eigenstates 
that, though having small energy splitting, do not show 
strong site correlation of quanta as in the previous cases. 
From these results we expect that escaping to the resis- 
tive state will not prevent from the experimental obser- 
vation of QB excitations. 

TABLE L Escape times for metastable states in a single 
Josephson junction Tescape estimated by formula (|17p . and the 
tunneling time of the initial excitationj^'(O)) = |0, 1) between 
the two junctions Ttunnei estimated from energy splittings. 



I 


'^tunnel (^s) 


"Tescape (^s) 


20 


348 


42 


19 


1.8 X 10^ 


3.5 X 10^ 


18 


10.16 X 10^ 


503.2 X 10'' 


17 


2.3 X 10^ 


71.2 X 10" 


16 


366 


1.62 X lO'-' 



Another phenomenon that was not taken into account 
in our quantum model is decoherence. To be able to ob- 
serve tunneling between the junctions the coherence time 
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has to be longer than the shortest tunnehng time between 
the junctions, which is on the order of 1 ns in the cases 
shown in Fig|4]-a and e. In the experiment shown in [H^l 
using a few levels per junction they obtained a coherence 
time on the same order. However, in the experiment in 
the coherence time was about 25 ns, and more recently 
in [3] the coherence time was approximately 80 ns. We 
expect that further improvements in experiments Q will 
give us even longer coherence times. 

Note that the above coherence times are shorter than 
the tunneling times of QB excitations (see table |T|, hence 
decoherence is an effect that can not be ignored if one 
wants to do a more realistic quantum description of the 
system. When exciting a JJ to high-energy states re- 
laxation (over dephasing) is usually the main source of 
decoherence. We can make a crude estimation of the cor- 
responding relaxation time Ti by using Ti ~ hQ/ei {Q 
is the quality factor of the junctions), which holds for a 
harmonic potential (tT], [T^I- For I — 19, 18 and 17, ei/h 
is around 150 GHz (see Fig[3]-b). For the JJs used in Q, 
Q is between 500 and 1000, which leads to a relaxation 
time between 3 ns and 6 ns. It is much smaller than 
the tunneling time of the QB excitations, therefore one 
would expect to see instead of tunneling, a freezing of the 
QBs on one of the junctions before they decohere due to 
relaxation. 

One could obtain more feasible results by increasing 
the bias current in such a way that there are less energy 
levels in the junctions. With this, exciting a QB would 
need less energy, and the relaxation time becomes longer. 
The tunneling time of that QB excitation is shorter, and 
might be even shorter than the relaxation time, allowing 
one to observe tunneling before relaxation. This possi- 
bility, and the inclusion of decoherence in our model, are 
issues that will be addressed in a future work. 



VII. CONCLUSIONS 

We have studied the classical and quantum dynamics 
of high-energy localized excitations in a system of two 
capacitively coupled JJs. In the classical case the equa- 
tions of motion admit time periodic localized excitations 
(discrete breathers) which can be numerically computed. 
For the quantum case we showed that excitation of one 
of the junctions to a high level, leaving the another junc- 
tion in the ground state, may strongly overlap with QBs 
(tunneling-pairs eigenstates) that live in the central part 
of the energy spectrum and localize energy on one of the 
junctions for a long (tunneling) time. This result would 
not qualitatively change if we excite a (coherent or in- 
coherent) superposition of several product basis states 
instead of only one. By using the density function for 
asymmetric superpositions of QB states one can realize 
how many quanta can be excited at each junction in order 
to excite QB states. 

In addition to what was described above, the sys- 
tem showed other interesting properties: We found that 



the system nearly conserves the total number of quanta, 
which comes from the fact that the coupling between 
the junctions couples just slightly eigenstates components 
with different total number of quanta. This opens the 
possibility to study larger systems without too big com- 
putational cost. When computing the fluctuation in the 
total number of quanta for each eigenstate, QB states 
show the smallest fluctuations. We showed that entan- 
glement, which reflects how many basis states have sig- 
nificant weight in an eigenstate, increases with energy in 
most of the eigenstates, becoming maximum at the center 
of the spectrum and then decreases. QB states from the 
same energy region are less entangled. This is because 
a QB state mainly consists of a symmetric or antisym- 
metric combination of two product basis states localizing 
many quanta on one of the junctions. 

With the available techniques for manipulating 
Josephson-junction qubits the experimental observation 
of QB excitations is possible. Escaping to the resistive 
state of the junctions (which together with decoherence 
was not taken into account in our quantum model) would 
not prevent us from doing that, and we expect that im- 
provements on preparation (higher quality factors) and 
isolation techniques of JJs make it possible to reach long 
enough coherence times in order to observe the phenom- 
ena we described in this work. By changing the param- 
eters of the system (bias current and coupling strength) 
one could vary the energy, and hence the tunneling time 
of a QB excitation with respect to the coherence time, 
in such a way that it becomes larger than that tunneling 
time. 
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APPENDIX: MEASURES OF ENTANGLEMENT 

Let A" ^ be a real eigenstate written in a ba- 
sis of product states {|ni,n2)} = {\ni) (g) |n2)} with 
ni,n2 = I, ■ ■ ■ N . It is a X iV matrix with elements 
Xni.n2 — ("ii^2|x)- We define the geometric measure of 
entanglement of this eigenstate by the following quantity: 

N 

^ = X! (Xnun2 - /r^l5r^2)^ (A.l) 

Where the vector functions f = (/i, . . . , /„^, . . .)* and 
g = (gi, . . . , (7„2J • ■ •)* such that A is minimum. The 
quantity A measures how far the eigenstate is from be- 
ing a product of functions depending respectively on the 
numbers ni and n2- 
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The minimization of A leads to the formula [6^. [65l. [6^ : 

A = llA-^f -A,„a., (A.2) 

where X'^ is a. N x N matrix with elements Xni n-i i "^max 
is the maximum eigenvalue of the N x N reduced density 
matrix 

A = X''iX''Y, (A.3) 

Another standard measure of entanglement in the 
eigenstates is the von Neumann entropy, which is used 
in information theory [69| : 

= -Tr {p, log, (A.4) 

= -^Afclog2(Afc), (A.5) 
fc 

where log2 refers to the logarithm taken in base 2. pi is 
the reduced density operator of either of the subsystems: 

Pi^Tr^ip), (A.6) 

where Tr, is the partial trace over the subsystem 2. {Xk} 
is the set of eigenvalues of the reduced density operator 
Pi- 

For the system of coupled JJs an eigenstate has the 
form: 

N 

\X) = X! X»l,n2l"l:"2)- (A.7) 
ni ,n2 

Then the density operator is 

P = IxXxl (A.8) 

N N 

= J2 Xni,n2Xn[,n'Jni,n2){n[,n'2\, (A.9) 



hence the reduced density operator is 



pi = Tr2ip) (A.IO) 

= E E^-.-^ni.«Jl"l>KI- (A.ll) 



The matrix elements of this reduced operator are 



N 

{n\pi\m) = Xl^"'"2^™>»2' (^-12) 

= ^n,m) (A. 13) 



where An_m is a matrix element of A defined in eq. (jA.Sp . 
To compute the von Neumann entropy one diagonalizes 
this matrix and uses eq. (|A.5p . 

Despite the fact that we found similarity in the varia- 
tion of the two measures when studying QB eigenstates, 
it is interesting to note, that monotonicity does not hold 
in general, i.e. if one measure is telling that a given state 
is stronger entangled than another one, that property 
may be reversed when using the other measure. The ge- 
ometric measure is an unambiguous number of the short- 
est distance from a given state to the subspace of prod- 
uct states, it allows to reconstruct the optimum product 
state, and it has a clear relation to the Hartree approxi- 
mation [g^l- For these reasons we presented the numeri- 
cal results using the geometric measure, rather than the 
entropy measure. 
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